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Viscoelastic models of ice-shelf flexure and ice-stream velocity perturbations are combined into a single 
efficient flowline model to study tidal forcing of grounded ice. The magnitude and timing of ice- 
stream response to tidally driven changes in hydrostatic pressure and/or basal drag are found to depend 
significantly on bed rheology, with only a perfectly plastic bed allowing instantaneous velocity response 
at the grounding line. The model can reasonably reproduce GPS observations near the grounding zone of 
Bindschadler Ice Stream (formerly Ice Stream D) on semidiurnal time scales; however, other forcings such 
as tidally driven ice-shelf slope transverse to the flowline and flexurally driven till deformation must also 
be considered if diurnal motion is to be matched. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

The influence of ocean tides on ice shelves and their tribu- 
tary ice streams provides opportunities for investigation of fun- 
damental processes controlling flow at ice-sheet margins. Early 
observational (Robin, 1958) and theoretical (Holdsworth, 1969, 
1977, 1981) studies of vertical ice-shelf flexure led to more com- 
plex modeling studies (Schmeltz et al., 2002; Reeh et al., 2003; 
Sergienko, 2010; Sayag and Worster, 2011, 2013), and ultimately 
to the use of remote sensing to detect the grounding zone over 
which ice transitions between grounded and freely floating states 
(Scambos et al, 2007; Fricker et al., 2009; Brunt et al., 2010, 
2011; Bindschadler et al., 2011; Rignot et al., 2011). Meanwhile, 
GPS observations demonstrated that the flow of several Siple Coast 
ice streams into the Ross Ice Shelf (Anandakrishnan et al., 2003; 
Bindschadler et al., 2003; Winberry et al., 2009, 2011) and the 
flow of Rutford Ice Stream into the Ronne Ice Shelf (Gudmundsson, 
2006, 2007; Murray et al., 2007; King et al., 2010) are modu- 
lated by the tides. Modeling studies of these velocity perturbations 
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suggest that the ice-stream beds can have markedly different char- 
acteristics, with Gudmundsson (2011) finding a basal flow law 
for Rutford Ice Stream consistent with hard-bed sliding or low- 
stress-exponent till deformation while Walker et al. (2012) find 
an effectively plastic basal flow law for Bindschadler Ice Stream 
(on the Siple Coast) consistent with sliding over weakly velocity- 
strengthening till (Tulaczyk, 2006). 

The precise mechanisms by which tidal motion causes modu- 
lation of ice-stream flow are not yet well known, but it is likely 
that improved understanding of these mechanisms will provide 
insight to grounding-zone processes that affect ice-stream dynam- 
ics. The two most obvious processes are the change in hydro- 
static pressure with tidal height and transient changes in ground- 
ing due to flexure, although modeling of the latter in an elas- 
tic, non-hydrostatic manner remains an active area of research 
(e.g., Sayag and Worster, 2011, 2013). Other potentially signif- 
icant processes include changes in ice-shelf slope due to spa- 
tially nonuniform tides (Brunt, 2008), flexurally induced variation 
in subglacial pressure affecting water flow (Walker et al., 2013), 
and formation of “sticky spots” due to repeated vertical com- 
paction of till (Christianson et al., 2013). Anandakrishnan and Al- 
ley (1997) were able to explain tidally forced basal seismicity of 
ice stream C (now Kamb Ice Stream) using a simplified model 
of an elastic ice stream with no explicit shelf overlying a viscous 
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bed. More recently, Gudmundsson (2011) used a viscoelastic, full- 
Stokes model to show that the fortnightly modulation of Rutford 
Ice Stream can be explained by interaction between semidiurnal 
tidal constituents provided that the basal sliding law is nonlinear. 

In this study, we combine and expand earlier models of ice- 
shelf flexure (Walker et al., 2013) and ice-stream response to stress 
perturbations (Walker et al., 2012) to analyze how tidal forcing 
leads to modulation of ice-stream flow, with an emphasis on ve- 
locity perturbations in the grounding zone. Beginning with exper- 
iments on idealized domains, we demonstrate the importance of 
basal rheology in determining the magnitude and timing of the 
response to a given forcing. The model is then applied to Bind- 
schadler Ice Stream, resulting in broad agreement with the obser- 
vations of Anandakrishnan et al. (2003). 

2. Model derivation 

The vertical displacement w of the ice shelf/stream (at time 
t (s) at a horizontal position x (m) along a chosen flowline) is cal- 
culated from the model of Walker et al. (2013), 


In Walker et al. (2012), which was concerned only with grounded 
ice, it was assumed that neither the basal pressure Pt, d x b nor the 
basal drag coefficient fi 2 varied in time; however, this assumption 
is no longer valid when vertical flexure of the ice shelf and/or 
tidally-driven changes in grounding are considered. For the basal 
drag term, we now have 

9f — &t (u + u ) 171 ) — — (u + u) m 1 9f U + 9f + u ) 171 , 

(7) 

which adds a term to the earlier model. As for the basal pressure, 
the elevation of the base of an ice shelf subject to tidal flexure is 
z = b(x, t) =Zb(x) + w(x, t), and its depth is r)(x , t) — b = r\ — Zb — w, 
where r] is the tidal height, w is the vertical flexure of the shelf 
(calculated from (1)), and z & is the initial basal elevation when r] = 
w = 0. We define z = 0 at sea level (i.e., the sea surface elevation 
when r] = 0). Assuming hydrostatic pressure at the ice-shelf base 
and using the preceding expressions for basal depth and elevation 
leads to 


Ek E 

d t (kw + 9* (D 9* w)) H =-w = d t q H (1) 

v xV x n 2v(l — A, 2 ) 4 2v(l — A. 2 ) 4 

which represents the ice as a viscoelastic beam with Young’s mod- 
ulus £, viscosity v, and Poisson’s ratio X resting (where grounded) 
on an elastic foundation with spring constant k. The flexural rigid- 
ity of the ice is given by 


D = 


Eh 3 

12(1 - A. 2 )’ 


( 2 ) 


where the thickness h can vary along flow. Hydrostatic pressure 
provides an applied load, 


q = Pwg0 7-w), (3) 

wherever the vertical displacement of the shelf differs from the 
tidal height r ]. The vertical displacement calculated by the flex- 
ural model comprises part of the forcing of the flowline model 
described below. 

To derive the perturbed flowline model (following Walker et al. 
(2012)), we first remove the assumption of viscous rheology from 
the Dupont and Alley (2005) momentum equation to obtain 


P b d x b = p w g(t]-z b -w)d x (z b + w), (8) 

so that 


3 t(P b dxb) = p w g{d t (r) - w)d x (z b + w) + (r]-z b - w)d x 3 t w j. 

(9) 


As in the earlier model, this term remains zero for grounded ice. 

We continue to take ice to be a Maxwell material, the simplest 
viscoelastic rheology that displays both instantaneous elastic re- 
sponse and long-time viscous behavior. This can be thought of as 
a spring and dashpot in series, so that the same stress acts across 
each element and the total strain e is the sum of the elastic and 
viscous strains, leading to the differential equation 


E E 

<7 = Ee < 7 = Ed x u cr, 

2v 2v 


( 10 ) 


(e.g., Turcotte and Schubert, 2002) where E is Young’s modulus 
and the viscosity v, given by Glen’s flow law, is assumed to change 
sufficiently slowly that the viscosity of the background state can be 
used in calculating the perturbations. The final perturbed momen- 
tum equation is obtained by substituting (7), (9), and (10) into (6): 


d x (2ha-^pgh 2 J = P b d x b + r b (u), (4) 

where u is velocity and a is along-flow deviatoric stress. The 
along-flow component of basal pressure Pt, d x b (which results from 
depth-integration) is the product of the pressure Pt, (cryostatic or 
hydrostatic) and slope at the bottom surface of the ice, while the 
basal drag under grounded ice is given by the power-law bed rhe- 
ology Tb(u ) = /3 z um, which can range from linear viscous (m = 1) 
to perfectly plastic (m -> oo) with the drag coefficient fi 2 (x, t ) typ- 
ically determined empirically. If we then consider small changes in 
velocity and stress, the perturbed momentum equation is 

dx^2h(a + or) — ^pgh 2 ^j = P b d x b + z b (u + u), (5) 

where a, u are the perturbation stress and velocity and a, u are 
the background state (derived from observations and/or standard 
viscous ice-flow modeling as in Walker et al. (2012)). Assuming 
negligible change in ice thickness or background state over the 
relatively short (on the order of days) timescales of interest, dif- 
ferentiating with respect to time gives 

d x (2h d t cr) = 3 t (P b d x b ) + d t (r b (u + «)). (6) 


9v 2hE 


«-*}) 


= Pwg{d t (t, - W) 3 x (rj — z b — w) 


+ (r,-z b -w)d x d t (tj-w)} 

+ — (u + U)m _1 3 t U + 3 t(/S 2 )(U + u)m , 
m y ’ 


HD 


where advective terms from the total time derivative in (10) are 
negligible. 

We discretize and solve (1) and (11) by the finite element 
method with fully implicit time differencing, using Hermite cubic 
and linear elements, respectively. For all experiments in this study, 
model resolution is 200 m in x and 600 s in t. 


3. Idealized experiments 

Our idealized flowline is 200 km long, consisting of a 150 km 
ice stream and a 50 km ice shelf. (As will be seen below, a longer 
shelf is not necessary for the simple experiments of this section.) 
The ice stream has a uniform thickness of 500 m, with the shelf 
tapering linearly to 400 m at the ice front. Experiments in which 
the shelf instead tapers following a concave ice-tongue profile 
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Fig. 1 . Perturbation velocity u at grounding line for 50 cm diurnal tide with various ice-stream bed exponents m. Negative of the rate of tidal height change (— 3 f rj) shown 
for comparison of timing. 


(Holland et al, 2008), leading to lower rigidity, produce negligible 
differences in vertical flexure (~1 mm for 50 cm tide). The ice- 
stream bed has a spring constant k of 10 7 Pam -1 , roughly in the 
middle of the reasonable range derived in Walker et al. (2013) us- 
ing C. Marone’s unpublished data from experiments in Rathbun et 
al. (2008). The ice has Young’s modulus E either 3.0 or 4.8 GPa 
and Poisson’s ratio X of 0.4 (Anandakrishnan and Alley, 1997; 
Gudmundsson, 2011). For most of our experiments in this sec- 
tion, the ice will have viscosity v = 10 20 Pas, a deliberately high 
value chosen to make elastic behavior dominant. The steady-state 
background velocity u increases linearly from 350 ma -1 at the up- 
stream end of the ice stream to 500 ma -1 at the grounding line, 
and the drag coefficient ff 2 is set for each exponent m so that 
the steady-state basal drag (calculated from the power-law bed 
rheology) decreases linearly from 50 kPa upstream to 20 kPa at 
the grounding line. We note that the latter value is likely on the 
high side for a grounding zone, making our idealized configuration 
rather sensitive to grounding-line motion. 

3.1. Tidal forcing with fixed grounding line 

In order to consider only the effects of hydrostatic pressure (i.e., 
the term given by (9)), we run these experiments with w = 0 
imposed at a fixed grounding line. The remaining boundary con- 
ditions are w = d x w = 0 at the upstream end and w = rj(t) at 

the downstream end; the initial conditions are w = r] = 0 every- 

where. For a diurnal tide with 50 cm amplitude, shelves with 
E = 3.0 and 4.8 GPa both have maximum vertical displacement 

w = 52.2 cm at high tide, with the stiffer (4.8 GPa) shelf having 

a somewhat longer flexure zone (11.0 vs. 9.8 km until |w — rj\ < 
1 mm). We note that flexure exceeding the boundary condition on 
vertical displacement (here, the tidal height) is common in beam- 
bending problems; cf. formation of forebulges during flexure of 
lithospheric plates. 

It might be expected that spreading the hydrostatic pressure 
change due to the difference between vertical flexure and tidal 
height over a wider grounding zone, thus decreasing the pressure 


gradient, would produce less forcing on the grounded ice. How- 
ever, the forcing is dominated by the hydrostatic pressure change 
at the grounding line, and thus the flexural parameters and the re- 
sulting shape of the flexure zone have negligible effect on the force 
felt by the grounded ice when the grounding line remains station- 
ary. We note that these parameters do become important when 
tidally driven grounding-line motion is calculated (e.g., Sayag and 
Worster, 2011, 2013), and that Young’s modulus of the grounded 
ice is a leading control on the magnitude of the velocity perturba- 
tion resulting from a given tidal forcing. Also, ice-shelf slope due 
to spatially varying tides can become significant when considering 
more realistic scenarios; this effect will be discussed later, when 
the model is applied to Bindschadler Ice Stream. 

The velocity perturbation at the grounding line resulting from 
a 50 cm diurnal tide is shown in Fig. 1 for E = 4.8 GPa and 
bed rheologies ranging from linear (m = 1) to perfectly plastic 
(m oo ). We note that because the ice is in a high-viscosity, 
elastic-dominated regime, the velocity is linear in Young’s mod- 
ulus, so the E = 3.0 GPa response is 1.6 times larger though iden- 
tical in timing. 

The case of a perfectly plastic bed, while unrealistic, is a useful 
starting point for understanding the model’s behavior. As m ->> oo 
(effectively, when m > 10 5 ), the second term on the RHS of (11) 
is eliminated, taking with it any effect of u on basal drag. (The 
third RHS term is also zero, as in this case we are not considering 
variation of the drag coefficient fi 2 due to grounding-line motion 
or basal hydrology.) With the perturbation stress term a /2v neg- 
ligible due to high viscosity, (11) reduces to a Poisson equation 
in u with time-dependent forcing. The solution is a perturbation 
velocity profile (linear in this case, as h is constant for the ice 
stream) in which all points respond instantly to the forcing. Thus, 
u at the grounding line is 180° out of phase with the time deriva- 
tive of the tide, 3 t r] (i.e., perfectly in phase with This result 

is to be expected, because for the grounded ice this simplified 
case is equivalent to applying a periodic load to an elastic body, 
for which displacement is linear in stress and velocity linear in 
stress rate. We note that this dependence of perturbation velocity 
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Fig. 2. Perturbation velocity u at grounding line for 25 cm semidiurnal tide with various ice-stream bed exponents m. Negative of the rate of tidal height change (—3 t r}) 
shown for comparison of timing. 


on the rate of applied stress strongly affects the magnitude as well 
as the timing of ice-stream response, with higher-frequency tidal 
constituents having a stronger effect than their amplitudes alone 
would suggest. For example, a semidiurnal tide (such as S2, 12.00 h 
period) will produce a stress rate and resulting velocity perturba- 
tion of comparable magnitude to a diurnal tide (such as Kl, 23.93 h 
period) with twice the amplitude (Fig. 2). 

As m decreases and the ice-stream bed becomes less than per- 
fectly plastic, the basal drag-velocity feedback takes effect, and (11) 
becomes a nonlinear diffusion equation with time-dependent forc- 
ing. Even at the grounding line, the ice stream no longer responds 
instantly to the rate of tidal height change. For a linear viscous 
bed (m = 1), u lags —d t rj by 3 h for a diurnal tide. This lag scales 
with the tidal period (e.g., 1.5 h for a semidiurnal tide), repre- 
senting a phase delay of 45°. (For convenience, our diurnal tide 
in these experiments has a period of exactly 24 h.) As —d t rj is 
always 90° behind r) for a sinusoidal tide, the lag between tidal 
extrema and the corresponding velocity extrema is 3/8 of the tidal 
period. The lag changes relatively little as m increases, shortening 
by only ~10 min for a diurnal tide when m = 100. Thus, for most 
realistic values of m, including the range (m e [8,40]) of “effec- 
tively plastic” values that Walker et al. (2012) found reasonable for 
Bindschadler Ice Stream, the bed rheology affects the magnitude 
but not the timing of the response of an elastic ice stream to tidal 
forcing. 

3.11 Effects of viscosity 

From the Maxwell rheology (10), the viscous and elastic contri- 
butions to the strain rate d x u are equal when a /a = 2 v/E, that is, 
when the time scale of the forcing equals the Maxwell relaxation 
time of the material. Forcing on much shorter time scales produces 
elastic behavior, and forcing on much longer time scales produces 
viscous behavior. With E = 4.8 GPa and v = 10 20 Pas, the relax- 
ation time is over 1300 years, confirming that the ice stream in 
our experiments thus far is in the elastic regime. In order to reduce 
the relaxation time to a day, the depth-averaged effective viscosity 
must decrease to approximately 2 x 10 14 Pas, a low but reasonable 
value for an ice shelf (Brunt, 2008). 


We re-run the diurnal tide experiment ( E = 4.8 GPa) with a 
perfectly plastic bed, so that any delay in ice-stream response is 
due solely to the effect of viscosity. In order to apply consistent 
forcing, the original (elastic) flexure results are used for all experi- 
ments, and only the flowline model is re-run. Results are shown 
in Fig. 3. For v ^ 10 16 Pas, the results are only negligibly dif- 
ferent from the elastic case, and for v = 10 15 Pas only a slight 
delay (^10 min) is noticeable. A larger delay (~20 min) occurs for 
v = 10 14 - 5 Pas, though the peak velocity due to the viscous contri- 
bution remains low enough (~16 ma -1 ) that the overall velocity 
does not noticeably increase. (Note that the strain rate, and thus 
the velocity, of a Maxwell material is the sum of viscous and elastic 
components, so that we can find the viscous contribution to each 
run by subtracting our earlier elastic results.) As the viscosity drops 
further, to v = 10 14 Pas, viscous effects become more readily ap- 
parent. The peak velocity due to the viscous contribution increases 
to ~ 51 ma -1 , increasing the overall peak velocity by 6.0 ma -1 , 
while the delay increases to 70 min. For v = 10 135 Pas, the vis- 
cous contribution (~157 ma -1 ) becomes comparable to the elastic 
contribution, increasing the overall peak velocity to 221 ma -1 and 
causing a lag of 3.0 h behind the elastic case. When v is reduced 
to 10 12 Pas, viscous effects dominate, leading to a peak veloc- 
ity over 2100 ma -1 offset by 6.0 h from the elastic case (not 
shown due to vertical scale). This timing indicates that the ice- 
stream response is synchronous with the applied stress, consistent 
with purely viscous behavior. However, the velocity is unrealisti- 
cally high when compared with observations (e.g., over an order 
of magnitude greater than seen by Anandakrishnan et al. (2003) 
at Bindschadler Ice Stream), even considering the also unrealistic 
assumption of a perfectly plastic bed. For most reasonable values 
of viscosity, we expect primarily elastic ice-stream response, with 
the possibility of a relatively small but noticeable viscous contri- 
bution. 

3.2. Grounding-line motion without tides 

To examine the effect of the basal drag terms (7) in isolation, 
we run experiments in which there is no tide, but grounding- 
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Fig. 3. Perturbation velocity u at grounding line for diurnal tide with perfectly plastic ice-stream bed and various ice viscosities v. Velocity for v = 10 20 Pas is identical to 
the m = 10 5 case in Fig. 1. 



Fig. 4. Perturbation velocity u at grounding line for ±100 m diurnal grounding-line motion with various ice-stream bed exponents m. Negative of the rate of change of 
grounding-line position (— d t Ax), i.e. rate of grounding-line retreat, shown for comparison of timing. Note 0.5 day shift of horizontal axis relative to Fig. 1. 


line motion of ±100 m (half an element) is imposed with the 
timing of a diurnal tidal cycle. (Although the present model can 
be parameterized to estimate grounding-line motion, we do not 
present results here because this process is better handled as a 
free-boundary problem (Sayag and Worster, 2011, 2013).) Thus, 
velocity perturbations are driven by the d t fi 2 term of (7), with 
the hydrostatic pressure term (9) absent. This forcing is of sim- 
ilar magnitude to the 50 cm tide of the earlier experiments 
as the change in horizontally-integrated basal drag is r^Ax = 
20 kPa • 100 m = 2 x 10 6 Pam, while for tide experiments, the 


change in depth-integrated hydrostatic pressure is approximately 
-p w grjZb = 1028 kg ITT 3 • 9.8 ms -2 • 0.5 m • -446 m ~ 2.25 x 
10 6 Pam. Results, shown in Fig. 4, are very similar to those of 
the diurnal tide experiments, with velocities proportional to the 
relative forcing (ratio ~0.89) and identical lags for a given bed ex- 
ponent m. Thus, the dependence of the ice-stream response on bed 
rheology is identical for the two forcings. It should be noted, how- 
ever, that the forcings themselves are 12 h (180°) out of phase and 
therefore act opposite each other. 
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Fig. 5. Bindschadler Ice Stream, showing locations of GPS stations (KO, K40, and K80) used by Anandakrishnan et al. (2003). Ross Ice Shelf grounding zone (Rignot et al., 
2011) in black, modeled flowline in white. Image from MODIS Mosaic of Antarctica (Haran et al., 2006). 


4. Application to Bindschadler Ice Stream 

We now apply the model to analyze data collected by Anan- 
dakrishnan et al. (2003) on Bindschadler Ice Stream (formerly Ice 
Stream D), which flows into the Ross Ice Shelf. Kinematic GPS mea- 
surements were taken at stations just downstream of the ground- 
ing zone (K0), 40 km inland (K40), and 80 km inland (K80) that 
lie near a best- fit flowline (Fig. 5). Analysis of the upstream prop- 
agation of tidally driven velocity perturbations by Walker et al. 
(2012) showed that the bed of this ice stream is effectively plastic 
(Rathbun et al., 2008), with the exponent m most likely between 8 
and 40. In the present study, we will focus on determining how the 
observed velocity perturbation in the grounding zone arises from 
tidal forcing, a question left open by the earlier work. 

The grounded section of our flowline is identical to that used by 
Walker et al. (2012), including flow parameters derived in 

that study by fitting surface velocity observations. The full flowline 
eventually makes a sharp northward turn, reaching the ice front 
just to the west of Roosevelt Island; because our one-dimensional 
model does not consider lateral drag, we include only the float- 
ing section of our flowline that extends essentially straight WNW 
from the grounding zone. We thus have a 438 km long flowline 
consisting of a 288 km ice stream and a 150 km ice shelf. Numer- 
ical resolution remains the same as for our idealized experiments 
(200 m in x, 10 min in t). 

Tidal forcing for our model is obtained from the Ross_Inv_2002 
regional tidal model (Padman et al., 2003), which has been opti- 
mized by assimilating gravimetry-derived tidal constituents on the 
Ross Ice Shelf. Only the I<1 (23.93 h), 01 (25.82 h), M2 (12.42 h), 
and S2 (12.00 h) constituents are included, because harmonic anal- 
ysis using the T_ TIDE (version 1.3 beta) program (Pawlowicz et al., 
2002) shows that these constituents dominate the velocity pertur- 
bations observed at station KO; each has a good signal-to-noise 
ratio, and our >11 day time series meets the Rayleigh criteria 
for separating constituents close in frequency. Given the strong 


vertical motion observed by Anandakrishnan et al. (2003), sta- 
tion K0 was certainly floating; however, given its proximity to the 
grounding zone and its smaller amplitude relative to modeled tidal 
heights, it was likely not fully hydrostatic. We note that although 
K0 does not lie perfectly on our flowline (~14 km north), the dif- 
ference between modeled tides at this station and at the grounding 
line of our flowline is negligible (a few mm per constituent). The 
grounding zone of Bindschadler Ice Stream is expected to have 
minimal motion based on ICESat repeat-track analysis (Brunt et al., 
2010), which suggests that the upstream limit of flexure is con- 
strained to within 500 m in this area. Therefore, we initially focus 
on the tides as the sole forcing mechanism. 

For all of the runs shown here, we use E = 3.0 GPa and m = 20. 
While we do not have enough data to uniquely determine these 
parameters, this combination produces velocities comparable to 
observations and has a bed rheology consistent with Walker et al. 
(2012). 

We begin with experiments in which the tide along the en- 
tire ice-shelf section of the flowline is assumed to be the same as 
the tide near the grounding zone. The resulting u at the ground- 
ing line is shown in Figs. 6 and 7. As observed by Anandakrishnan 
et al. (2003), the perturbation velocity reaches its maxima on the 
falling tide and its minima on the rising tide. Due to varying re- 
sponse times for the tidal constituents, the lag between tidal and 
velocity extrema (using all four constituents together) has a mean 
of 6.59 h and standard deviation of 2.63 h. Observed lags at station 
K0 (calculated by cross-correlation of the along-flow displacement) 
were 6±1 h (Anandakrishnan et al., 2003), so the model results 
include the observed range but also some significantly different 
lags. To investigate this behavior, we consider results for individual 
tidal constituents. (While the velocity-dependent basal drag makes 
(11) nonlinear, comparison of results for constituents used individ- 
ually and in combination shows that the effect is weak enough to 
allow a term-by-term approach.) Model response to diurnal tides 
(9.02 h lag for Kl, 9.73 h for 01) is noticeably slower than observed 
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Fig. 6. Perturbation velocity u in grounding zone of Bindschadler Ice Stream for E = 3.0 GPa and bed exponent m = 20. Results are shown for spatially uniform tide and for 
tide varying along the 150 km floating section of the flowline. Tide rj shown for comparison of timing. For clarity, only every third point of u is plotted. 



Fig. 7. Perturbation velocity u in grounding zone of Bindschadler Ice Stream for E = 3.0 GPa and bed exponent m = 20. Results are shown for spatially uniform tide and for 
tide varying along the 150 km floating section of the flowline. Rate of tidal height change 3 t r) shown for comparison of timing. For clarity, only every third point of u is 
plotted. 


(7.05 ± 0.18 h, 8.39 zb 0.41 h respectively, calculated using the 95% 
confidence intervals from T_TIDE), but slightly quicker for semidi- 
urnal tides (4.71 h, 4.87 h for S2, M2 respectively vs. 5.58 zb 0.35 h 
observed for S2). We note that the lag for the M2 constituent 
cannot be reliably calculated from the observations because the 
apparent M2 amplitude of the vertical motion is smaller than the 
uncertainty of the GPS measurement. 

Due to the great size of the Ross Ice Shelf, however, variations 
in the phase of tidal height across its extent are potentially sig- 


nificant (e.g., Brunt, 2008). We next run experiments in which the 
variation of the tide along the ice shelf is considered, so that even 
the freely floating part of the shelf (where w = ij) will have some 
slope, and the second term of (9) can thus be nonzero over the 
entire shelf. (While a one-dimensional flowline model cannot fully 
capture this inherently two-dimensional situation, it is reasonable 
to expect the along-flow and transverse forcings to be mostly inde- 
pendent when far from lateral boundaries.) Although d x r] is rather 
small (tidal height difference <10 cm over 150 km), its effect 
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Fig. 8. Example of tidal height variation along the ice-shelf section of the flowline, shown at successive low and high tides. Tidal heights are calculated from the Ross_Inv_2002 
regional tidal model (Padman et al., 2003). 


over freely floating ice tends to be directed opposite the effect 
of the near-grounding-zone tide. That is, most of the shelf slopes 
downward along flow (leading to a downstream-directed force) 
when high tide is reached near the grounding zone (an upstream- 
directed force), and vice versa (Fig. 8). The result (Figs. 6 and 7) 
is a reduction in the magnitude of velocity extrema (typically 
8-12%) and additional lag (typically 10-30 min) compared to the 
uniform-tide run. Separate runs show that the lags increase for 
each constituent (9.33 h, 10.06 h, 5.25 h, 5.11 h for Kl, 01, M2, S2 
respectively), leading to a mean overall lag of 7.28 h with standard 
deviation 2.95 h, again including the observed range. 

The model thus does a reasonably good job of reproducing 
semidiurnal response, but the diurnal response in the grounding 
zone of the best-fit flowline is roughly 2 h later than observed 
at station I<0. As seen in our idealized experiments, reducing the 
lag by this much would require assuming that the ice-stream bed 
is highly plastic (m > 500) for diurnal forcing but has a drasti- 
cally lower flow exponent for semidiurnal forcing. A frequency- 
dependent flow exponent would be possible if till properties were 
depth-dependent, as the depth of interaction will be greater for 
diurnal than for semidiurnal forcing. However, deeper till and/or 
bedrock would be expected to be less plastic, opposite to the be- 
havior needed to explain the observations. Furthermore, the anal- 
ysis of upstream propagation of velocity perturbations by Walker 
et al. (2012) is most consistent with a single flow exponent much 
lower than would be required to match the lag in diurnal response. 
It appears more likely that at least one other forcing on a diurnal 
timescale is present at station K0, and that the timing of the ne- 
glected forcing(s) leads to an earlier total response. 

It is possible that the diurnal response is affected by the geom- 
etry of the grounding zone near station KO, which is significantly 
more complex than near the best-fit flowline (Fig. 5). In this re- 
gion, flowlines (corresponding to streak lines in the MODIS image) 
are not noticeably influenced by the irregularities of the grounding 
zone, suggesting that the bed is not locally strong enough to affect 
the direction of the mean flow. The streak line through KO does 
not run seaward along the inlet; instead, it regrounds downstream 


of the station, and passes through areas of significant ephemeral 
grounding both upstream and downstream. 

Anandakrishnan et al. (2003) observed significant tidally driven 
motion of station KO transverse to the mean flow of the ice stream, 
speculating that this cross-line perturbation was due to the same 
shelf tilt that we consider only in the along-flow direction. We per- 
form tidal analysis of these GPS velocities at KO using T_TIDE. 
While GPS-related systematic errors associated with the Kl and 
S2 constituents probably exist, we assume that they are generally 
small (<3 mm) based on King et al. (2008, 2011). Our results in- 
dicate a complicated two-dimensional velocity vector composed of 
ellipses associated with the tidal constituents (Fig. 9), with a strong 
tendency for the velocity perturbation to be aligned with the in- 
let or the upstream area of ephemeral grounding rather than with 
the flowline. The complexity of this forcing, which may involve the 
entire shelf, could in itself lead to results that cannot be predicted 
by a one-dimensional model. Furthermore, this combination of off- 
axis forcing and ephemeral grounding should lead to constantly 
shifting patterns of basal drag along the direction of the perturba- 
tion and thus affect the timing of the response at station KO. Yet 
another possibility is that flexure-driven pressure changes at the 
bed (Walker et al., 2013) cause subglacial water flow and/or till 
deformation that could also produce fluctuations in basal drag. 

For a given tidal constituent, any process affecting basal drag 
would likely result in time dependence of the lag and/or of the 
amplitude of the velocity perturbation. However, harmonic analy- 
sis programs like T_TIDE seek constant-amplitude sinusoidal sig- 
nals at known tidal frequencies, so that any modulation of the 
grounding-line response (in frequency or amplitude) present in 
the data would be averaged out in the fitting process. Detection 
of grounding-zone basal drag effects may then depend on other 
signal-processing techniques like bandpass filtering; unfortunately, 
the current time series does not have sufficient length or reso- 
lution to produce unambiguous results. Further investigation will 
require a two-dimensional model, along with more detailed obser- 
vations of the relative timing of tides and ephemeral grounding 
throughout the grounding zone. 
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Fig. 9. (a, b, c) Observed velocity at station KO for Kl, 01, and S2 tidal constituents, respectively. One tidal cycle (low to low) shown for each constituent, (d) Observed velocity 
at station KO for Kl, 01, and S2 tidal constituents combined. Days 342 to 352 of year 2000 shown. Velocity vector follows path in a generally clockwise direction. 
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5. Conclusions 

Understanding the mechanisms by which vertical flexure of ice 
shelves induces horizontal motion of ice streams will likely lead 
to improved understanding of ice-stream dynamics. Our viscoelas- 
tic coupled flexure-flowline model shows that over the full range 
of bed exponents, the magnitude and timing of ice-stream re- 
sponse to tidal forcing depend strongly on the basal sliding law, 
although there is a range of realistic, effectively plastic bed expo- 
nents for which timing is nearly independent of rheology. Because 
we consider primarily diurnal and semidiurnal timescales, the elas- 
tic component of the response dominates, so that the ice stream 
responds to the rate of applied stress (with some delay as the slid- 
ing law departs from perfect plasticity). Applying this model to 
Bindschadler Ice Stream, we obtain reasonable agreement with ob- 
servations of highest velocity during the falling tide and lowest 
velocity during the rising tide. More precise matching of observed 
velocities will require both a higher-dimensional model and de- 
tailed observations of ephemeral grounding throughout the geo- 
metrically complex grounding zone. 
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